Tests of the Gravitational Inverse-Square Law below the Dark-Energy Length Scale 
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We conducted three torsion-balance experiments to test the gravitational inverse-square law at 
separations between 9.53 mm and 55 /im, probing distances less than the dark-energy length scale 
Ad ~ y^hc/pd ~ 85 /im. We find with 95% confidence that the inverse-square law holds (|a| < 1) 
down to a length scale A = 56 /^m and that an extra dimension must have a size _R < 44 fim. 

PACS numbers: 04.80.-y,95.36.+x,04.80.Cc,12.38.Qk 



Recent cosmological observations [H 0, S] have shown 
that 70% of all the mass and energy of the Universe 
is a mysterious "dark energy" with a density ~ 3.8 
keV/cm'^ and a repulsive gravitational effect. This dark- 
energy density corresponds to a distance Ad = y^c/pd ~ 
85 nm that may represent a fundamental length scale 
of gravity 0, 0|. Although quantum- mechanical vacuum 
energy should have a repulsive gravitational effect, the 
observed pd is between 10^*^ to 10^^" times smaller than 
the vacuum energy density computed according to the 
standard laws of quantum mechanics. Sundrum[6i] has 
suggested that this huge discrepancy (the "cosmologi- 
cal constant problem") could be resolved if the gravi- 
ton were a "fat" object with a size comparable to Ad 
that would prevent it from "seeing" the short-distance 
physics that dominates the vacuum energy. His scenario 
implies that the gravitational force would weaken for ob- 
jects separated by distances s < Ad- Dvali, Gabadaze 
and Senjanovicj?! argue that a similar weakening of grav- 
ity could occur if there are extra time dimensions. In 
their scenario, the standard model particles are localized 
in "our" time, while the gravitons propagate in the ex- 
tra time dimension(s) as well. Other scenarios predict 
the opposite behavior: the extra space dimensions of M- 
theory would cause the gravitational force to get stronger 
for s ^ i? where R is the size of the largest compactified 
dimension [s'l. These considerations, plus others involving 
new forces from the exchange of proposed scalar or vec- 
tor particles (§1 motivated the tests of the gravitational 
inverse-square law we report in this Letter. 

Our tests were made with a substantially upgraded 
version of the "missing mass" torsion-balance instrument 
used in our previous inverse-square-law tests [13, [HI ■ The 
instrument used in this work|12|. shown in Fig. [1] con- 
sisted of a torsion-pendulum detector suspended by a 
thin « 80-cm-long tungsten fiber above an attractor 
that was rotated with a uniform angular velocity w by a 
geared-down stepper motor. The detector's 42 test bod- 
ies were 4.767-mm-diameter cylindrical holes machined 
into a 0.997-mm-thick molybdenum detector ring. The 
hole centers were arrayed in two circles, each of which had 
21-fold azimuthal symmetry. The attractor had a simi- 
lar 21-fold azimuthal symmetry and consisted of a 0.997 



FIG. 1: Scale drawing of our detector and attractor. The 3 
small spheres near the top of the detector were used for a 
continuous gravitational calibration of the torque scale. Four 
rectangular plane mirrors below the spheres are part of the 
twist-monitoring system. The detector's electrical shield is 
not shown. 



mm thick molybdenum disc with 42 3.178-mm-diameter 
holes mounted atop a thicker tantalum disc containing 21 
6.352-mm-diameter holes. The gravitational interaction 
between the missing masses of the detector and attractor 
holes applied a torque on the detector that oscillated 21 
times for each revolution of the attractor, giving torques 
at 21a;, 42a;, 63a;, etc. that we measured by monitor- 
ing the pendulum twist with an autocoUimator system. 
The holes in the lower attractor ring were displaced az- 
imuthally by 360/42 degrees and were designed to nearly 
cancel the 21a; torque if the inverse-square law holds. On 
the other hand, an interaction that violated the inverse- 
square law, which we parameterize as a single Yukawa 

V{r) = -G^^^ [1 + a exp(~r/A)] , (1) 
r 

would not be appreciably canceled if A is less than 
the 1 mm thickness of the upper attractor disc. We 
minimized electromagnetic torques by coating the en- 
tire detector with gold and surrounding it by a gold- 
coated shield consisting of a tightly-stretched, 10 ^m- 
thick, beryllium-copper membrane between the detector 
and attractor plus a copper housing that had small holes 
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for the suspension fiber and the autocoUimator beam. 
The entire system was under a vacuum of ~ 10~® torr in 
a temperature-controlled and magnetically-shielded en- 
vironment. The noise in our torque measurements was 
generally close to the thermal value expected from the 
finite quality factor, Q « 3000, of our torsion oscillator, 
but increased noticeably for detector-membrane separa- 
tions below 100 fim (see Fig. [2]). A continuous, absolute 
calibration of the torque scale was provided by the grav- 
itational octupole interaction between 3 small spheres 
mounted on the detector and 3 larger spheres, mounted 
outside the vacuum vessel on a turntable that rotated the 
spheres about the fiber axis at a steady rate lUc , providing 
a calibration signal at 3a;c. 

The detector/attractor separation, ^ = {x, y, s), where 
X and y are the horizontal displacements between the 
centers of the detector and attractor and s is the vertical 
separation between the bottom of the detector ring and 
the top of the upper attractor disc, was determined using 
capacitance plus micrometer techniques for s, and grav- 
itational plus micrometer techniques for x and y. De- 
tector twist data were taken at x and y close to zero 
(except for off-center runs used to find the a;, y center) 
and separations 55 /um < s < 9.53 mm. The key pa- 
rameters, T]™'^ ± of our instrument (the masses 
removed in machining the holes, the hole radii, thick- 
nesses and locations, etc.) were determined precisely 
using an electronic balance or a coordinate-measuring 
machine (CMM). The CMM readings were corrected for 
surface-roughness by scanning the relevant surfaces with 
an atomic-force-microscope; the correction increased hole 
diameters and decreased plate thicknesses by 4.4 and 2.3 
/im, respectively. We measured the attractor angle, 0, 
by counting pulses to the stepper-motor and based (50 on 
the recorded scatter of the (p readings when a once-per- 
revolution index pulse occured. 

Our analysis strategy and techniques were described in 
detail in Ref. [ll|. Briefly, pendulum twist signals were 
digitally filtered to suppress the free-torsional oscillations 
and converted into torques by taking into account pen- 
dulum inertia and damping as well as signal-averaging 
and filtering. However, in this work, we made two im- 
provements to our analysis procedure. 

1) We used a more sophisticated 5-point torsion filter that 
also removed the effects of slow drifts in the equilibrium 
twist of the torsion pendulum, allowing us to eliminate 
the polynomial "drift terms" that were needed to fit the 
twist signals in Refs. 111- 

2) We mapped the linearity of our autocoUimator system 
after each run by stopping the attractor and calibration 
turntables and setting the pendulum into a free oscilla- 
tion that covered the same region of our photodetector 
as the preceeding data run. 

The resulting extracted torques were then decomposed 
into harmonic amplitudes at multiples of the attractor 
rotation frequency uj. 
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FIG. 2: Fourier transform of the raw twist signal in Experi- 
ment III taken at s = 67 (a detector-membrane separation 
of 46 ^m). The detector's free resonance occurs at 7.5tJ and 
the gravitational calibration is at Quj. The peaks at 21, 42 and 
Q'iuj probe the inverse square law. The smooth curve shows 
the thermal noise level. At this small separation the torque 
noise power retains the expected 1// form, but its amplitude 
exceeds the thermal value by about a factor of four. 



The observed harmonic torques, Nm{Cj) ± S^miCj), 
where j runs over all C values and to=21 or 42 (signals 
from higher m torques were significantly attenuated by 
pendulum inertia), were fitted by predicted torques, Nm 
that were functions of the key physical parameters of 
the apparatus, fj. We computed the Nm by integrat- 
ing Eg. [1] over our measured geometry as described in 
Ref. We first considered the effects of Newtonian 

gravity alone. We used gravity to determine indepen- 
dently the key physical parameters of the apparatus, ff, 
and required the two determinations to agree within er- 
rors. This was done by minimizing 



EE 

j m 



AA^m(O) 



-E 



cxp 
In 



Srj', 



cxp 



where the experimental error 



A^^(O) 



Ss 



dNrn 



(2) 



(3) 



accounted for the uncertainty (typically ±l/xm) in our 
measurement of the detector's vertical position. This fit- 
ting procedure yielded uncertainties that included statis- 
tical and most systematic effects, and took into account 
possible correlations between the uncertainties in the var- 
ious experimental parameters. Finally, we expanded the 
analysis to include a single Yukawa interaction. 

We conducted three separate experiments with our ap- 
paratus. These used the same detector ring and upper at- 
tractor disc, but different lower attractor discs. In Exper- 
iment 1 the lower-attractor disk thickness was 3.032 mm 
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FIG. 3: [Color online] Experiment I torques as functions of 
vertical separation s. The horizontal scale is expanded below 
s = 100 jj,m. The upper panel displays the 21w and 42tj 
torques as solid circles and squares, respectively. When not 
visible the errors are smaller than the size of the points. The 
smooth curves show the Newtonian fit to the combined data 
of all 3 experiments. The lower panel shows the 21tj residuals; 
the smooth curve shows the residual that would arise from an 
a = 1, A = 80 Yukawa interaction. 



and the Newtonian 21cj torque was over-canceled at all 
values of s. In Experiments II and III, this thickness 
was reduced by 0.140 mm, causing the Newtonian 21a; 
torque to be under-canceled for s < 100 /im. On the 
other hand, a short-range interaction with A << 1 mm 
would produce the same torques in all three experiments, 
but it must show opposite behaviors in Experiment I 
compared to Experiments II and III; decreasing, for ex- 
ample, the magnitude of the 21cl; signal in one case and 
increasing it in the other. This feature was useful in dis- 
criminating against possible systematic errors. After the 
two-plate attractor data were taken we accumulated ad- 
ditional data with just the upper or lower attractor plate. 

We used a 20 /xm-diameter suspension fiber in Ex- 
periment I, and set the attractor rotation frequency uj 
to ujo/28 and the calibration turntable frequency ujc to 
49uj/3, where coq « 12.6 mrads/s was the detector's free- 
oscillation frequency. After the Experiment I data were 
taken we discovered that the detector ring was slightly 
bowed so that its outer set of holes was slighter higher 
than the inner set. We accounted for the detector's de- 
formation by modeling the outer and inner sets of holes 
at different average heights, Az = 3.5 /im, above the at- 
tractor. (This was not necessary in Ref. pjj which used 
a thicker and more rigid detector.) Data and the New- 
tonian fit are shown in Fig. [31 A possible deviation at 
s < 80 /im led us to perform Experiments II and III. 
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FIG. 4: [Color online] Experiment II torques. Notation is the 
the same as for Fig. O but the vertical scales are expanded. 
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FIG. 5: [Color online] Experiment III torques. Notation is the 
same as in Fig. |4l except that diamonds (triangles) show the 
21u! torque from the upper (lower) attractor plate alone. The 
solid and dashed curves in the lower panel show the residuals 
expected from a = I, A = 80 /im and a = 10*, A = 10 /tm 
Yukawa interactions, respectively. Both are excluded by our 
results. 



Before beginning Experiment II we flattened the detec- 
tor ring to Az < 1.0 /fm, and switched to a 17 /xm fiber 
that gave ujq — 9.7 mrads/s. We reduced the 1// torque 
noise from internal losses in the suspension fiber by in- 
creasing the attractor and calibration turntable frequen- 
cies to a; = a'o/7.5 and uJc = 6a;, respectively. Results 
from this experiment are shown in Fig. 3) The small- 
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FIG. 6: [Color online] Constraints on Yukawa violations of 
the gravitational law. The shaded region is excluded at 
the 95% confidence level. Heavy lines labeled Eot-Wash 2006, 
Eot-Wash 2004, Irvine, Colorado and Stanford show experi- 
mental constraints from this work, Refs. [ij, [3], and 
01 [13 1 respectively. Lighter lines show various theoretical 
expectations summarized in Ref. [^. 



s anomaly is not evident. Hoping to reduce the excess 
noise at small s, we disassembled the apparatus before 
Experiment III and replaced the gold coatings on the 
detector and beryllium-copper membrane. In the pro- 
cess we inadvertantly increased the deformation of the 
detector ring to Az = 3.9 /im. We returned to a 20 ^m 
fiber, and operated at w = ujo/7.5 and cUc — 3uj. Data 
are shown in Fig. [5l The slightly different dependence 
of the 21uj torque on s, compared to Experiment II, is 
due to the detector curvature in Experiment III. In this 
case a marginally significant short-range anomaly corre- 
sponding to an increased attraction, similar to that seen 
in Experiment I, is again present. 

Non-gravitational backgrounds from magnetic, elec- 
trostatic and thermal effects were studied and found 
to be negligible using techniques described in Ref. 
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Casimir forces can only affect our results if they could 
penetrate the beryllium-copper membrane or cause the 
membrane to deflect so as to trace the hole pattern of the 
attractor. Private communications from Robert Jaffe, 
Paul Chesler, Anton Andreev and Laurence Yaffe have 
shown that our beryllium-copper membrane was thick 
enough to reduce direct Casimir forces between the at- 
tractor and detector to a negligible level. The observed 
1.6 kHz frequency of the lowest "drum-head" mode of the 
membrane was sufficiently high that m = 21 deforma- 
tions were insignificant. No corrections for backgrounds 



were necessary. 

A combined Newtonian fit to the data from all 3 ex- 
periments gave = 407 for v — A21 degrees of freedom. 
The best fit with an additional Yukawa interaction im- 
proved the by 3.5 for a = -0.0037, A = 2 mm. The 
combined data showed no evidence for a 2a effect at any 
A. Our resulting constraints on violations of the inverse- 
square law, shown in Fig.[6l improve on previous work by 
a factor of up to 100. In particular, at 95 % confidence, 
we find that any gravitational-strength [\a\ — 1) Yukawa 
interaction must have A < 56 /im. The results in Fig. [5] 
yield a model-independent upper limit on the size of a 
compact extra dimension. A single extra dimension with 
R ^ Smin would givc a signal corresponding to a Yukawa 
interaction with a = 8/3 and A = i? ^9], leading to a 
95%-confidence upper bound of i? < 44 /xm. For the two 
large extra-dimension scenario discussed in Ref. we 
require a 2a lower limit on unification mass M.^, > 3.2 
TeV/c^, where M* is defined in Ref. [lH. Constraints 
from the data in Figs. 13141 [5] on other possible forms of 
inverse-square-law violation will be submitted as a sepa- 
rate publication. 
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